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Abstract. - We extend the quantum geometric tensor from the state space to the operator level, 
and investigate its properties like the additivity for factorizable models and the splitting of two 
kinds contributions for the case of stationary reference states. This operator-quantum-geometric 
tensor (OQGT) is shown to reflect the sensitivity of unitary operations against perturbations of 
multi parameters. General results for the cases of time evolutions with given stationary reference 
states are obtained. By this approach, we get exact results for the rotated XY models, and show 
relations between the OQGT and quantum criticality. 
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Introduction. — The geometric structure on the 
manifold of quantum states (MQS) is an interesting sub- 
j^Jjpct of quantum physics, which was studied from different 
^^perspectives, e.g., the geometrization of quantum mcchan- 
C\) ips [U i the geometric phase [2H5] , the geometry of quantum 
CO evolution [5] , the quantum theory of gravity [THH] , quan- 
^^tium phase transitions [TUHH], etc. Generally speaking, 
00 on the MQS, there is a natural Riemann structure — the 
•O 1 quantum-geometric tensor (QGT), induced by the inner 
product in the Hubert space H 0H2]. The real part of 
. .the QGT provides the metric tensor, through which we 
can measure the distance between two states on the MQS. 
k>( The imaginary part of QGT gives a curvature 2-form, cor- 
<_i responding to a natural symplectic structure on the pro- 
spective Hilbert space V(H). With this curvature 2-form, 
the geometric phase can be interpreted as the area en- 
closed by the closed curve in V(Ji) [T3] . 

The corresponding metric given by the QGT is a kind of 
Fubini-Study metric related to the fidelity of two quantum 
states, which is one of the most popular physical quanti- 
ties in the quantum-information processing. Recently, it is 
proposed that the fidelity can characterize quantum phase 
transitions [14]. This approach has been achieved a lot 
of progress [TUl[TllH4H5n] . More recently, the fidelity ap- 
proach was extended to the operator level — the operator 
fidelity (OF) [3"Tll3"2]. The underlying idea behind the OF 
is mapping the unitary operators into the quantum states. 



It is proposed that the susceptibility (metric) of the OF 
can be used to study the stability of quantum evolution 
with respect to perturbations, indicate the quantum criti- 
cal points, quantify the environment induced decoherence 
and investigate quantum chaos |31H35j . 
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Since unitary operators can be mapped into quantum 
states, on the manifold of the unitary operators , we can 
define an operator-quantum-gcometric tensor (OQGT), 
whose real and imaginary parts can give a metric tensor 
and a curvature 2-form respectively. Like the QGT is nat- 
urally induced by the inner product in the Hilbert space 
%, the OQGT is naturally induced by the inner product 
of two operators, which is a generalization of the Hilbert- 
Schmidt product, with a reference state p [32]. Like the 
susceptibility of OF, The OQGT has the additivity for the 
factorizable model and the splitting of two kinds of contri- 
butions for the special cases that the reference state p is 
commutable with the operators. We will show this OQGT 
approach is useful to the sensitivity of unitary operations 
against multi-parameter perturbations. 



Riemannian structure on the MQS. The term 
MQS is pointed to the submanifold of the projective 
Hilbert space, which is defined as the sets of rays of the 
Hilbert space H. The inner product structure of the 
Hilbert space % naturally induce a Riemannian structure 
of the projective Hilbert space V{H) [T2]. We denote the 
ray of W by \tf>). The Fubini-Study metric indicating the 
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distance between \i[>) and \(j>) is denned as follows 

7 (ip, 4>) = arccos-F, (1) 

where F = \(ip\<p)\ is the fidelity defined through the inner 
product in W. The infinitesimal form of this metric, which 
denote the distance between two slightly different states 
\ip(X)) and \i/(X + dX)), can be written as 

ds 2 ~A(\-F 2 ) = Q py dX»dX\ (2) 

where 

Q^^V 1(1 -1^1)1 W (3) 

is the so-called quantum geometric tensor on V(H) [51IT2] . 
Here d p = d/dX^ with A M -s the parameters characteriz- 
ing the elements in T(H). Hereafter, we use the Einstein 
summation convention. Q pLI is a Hermitian matrix, whose 
real part is symmetric while the imaginary part is anti- 
symmetric. Hence, only the real part of quantum geomet- 
ric tensor has effects on Eq. Q due to the summation , 
i.e., (is = q^lu 

dX^dX", where 9nv — ReQ^j,. Meanwhile, 
the imaginary part gives the curvature 2-form 

(T = (ImQ M „) dX^ A dX", (4) 

where A is the exterior (wedge) product. The integration 
of (7 over an area S gives the geometric phase 7 g = — JJ S a 
of the cycle evolution along the boundary of the area S 
2,4 . Further, a can be written as the exterior differential 
of a 1-form Berry's connection (3 = —i(ijj\di]j), i.e., a = 
d/3 = -i(dip\ A \dip). 

Operator quantum geometric tensor. In the 

space of operators, one can define the following Hermitian 
inner product 

(X,Y) P :=Tr(XtY», (5) 

where p is a density operator taken as the reference state. 
Especially, we concern the unitary operators hereafter, 
which satisfies the normalized condition (U, U) p = 1. Ac- 
tually, this inner product structure can be interpreted as 
the conventional one for the two pure states \X) p and \Y) P 
in the extended Hilbert space |32j , through a map from the 
operator X to the pure state \X) P given as follows 

X^\X) p = X®l\*{p)), (6) 

where 1 is the identity operator in the ancillary 
Hilbert space "H anc and \^(p)) is a pure state satisfying 
Tr anc (\^ (p)) (^ (p)\) = p, i.e., it is a purification of p in 
the extended Hilbert space. Obviously, the concrete form 
of purification does not influence the inner product ([5]). 
For brevity, we omit the subscript label p of \X) p here- 
after as long as not confusing. 

Follow this spirit, the fidelity of two unitary operators 
is defined by 

F(U 1 ,U 2 ) = \Tr(ulU 2 p)\, (7) 



which is equivalent to the fidelity |(£/i|[/2)| of the two pure 
states \U\) and | C/2) - If U\ and U2 are slightly different 
and p is pure, F (1/1,1/2) is just the square root of the 
Loschmidt echo [36], which is used to describe the hyper- 
sensitivity of the time evolution to perturbations. It is 
remarkable that the minimization of F (1/1,1/2) with re- 
spect to a optimal p is used to characterize the statistical 
distinguishability of two unitary operations U\ and U2 |37j . 

Like the Riemannian structure on the MQS [12], we can 
naturally induce the Ricmann structure on the manifold of 
the unitary operators along with the given reference states. 
The main geometric object characterizes this structure is 
the quantum geometric tensor of the unitary operators 

Q pv = (d p U\d„U) - (d p U\U) (U\d v U) 

= (A fl A v ) p -(A fl ) p (A v ) p , (8) 

where A p = iU^d p U is the Hermitian generator of dis- 
placements in A M and (O) = Tr [pO] denotes the ex- 
pected value of the observable O on the state p. Hereafter, 
we call Q pl , operator-quantum-geometric tensor (OQGT). 
The real part g uv = Re Q pu is the Fubini-Study metric 
tensor of the unitary operator space with the reference 
state p, i.e., ds 2 = g uv dX^dX v , which measure the statisti- 
cal distance [551I3T?] between the two states U (A) ®l\^>(p)) 
and U(X + dX) ® l\^>(p)) in the extended Hilbert space. 
Meanwhile, due to the anti-commutative law of the ex- 
terior product, the imaginary part ImQ^^ = Im^^^j,) 
gives the curvature 2-form as follows 

a=-i(AAA) p , (9) 

where A = A p dX^ L = iU^dU is a differential 1-form. If 
the reference states p are independent on A^-s, a can be 
written as a = d/3, where (3 = —(A) p is the analog of 
Berry's connection. Because j3 can be written as (3 = 
-i(^>(p)\(W ®T)(dU ®T)\^>(p)) with *(p) the purification 
of p, it is, in fact, the Berry's connection in the extended 
Hilbert space. For a closed trajectory dS — {X(t) : t G 
[0,T)} with A(0) = X(T) in the parameter manifold, the 
geometric phase is given by 7 g = — <f> ds f3 = — JJ„ a, where 
5* is the area subtended by the closed trajectory dS. 

Next, we investigate the properties of the OQGT for 
the cases of factorizable models and of stationary reference 
states. These will be useful to the calculations in the XY 
model presented later. 

Additivity for factorizable models. We analyze the 
case in which the unitary operator U(X) and the refer- 
ence state p(X) have the same factorization structure as 
follows 

M M 
U(X) = ($U k (X), p(A) = (g)p fc (A). (10) 

This assumption is sensible when U(X) and p(X) are both 
generated by the same Hamiltonian. For instance, U = 
exp (-itH(X)) is the time evolution generated by H(X) 
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and p(A) is the Gibbs thermal state which is a mixture of 
the eigenstates of H(X). After derivative with respect to 
A M we obtain 

M l-l M 

d„U(\) = J2[(g)U k (X)}®dt l U l (\)®[(g) U k (\)}. (11) 

7=0 k=0 k=l + l 

Because Tr(A ® B) = Tr(A)Tr(5) and Tr(p z ) = 1, we get 



Note that [Aj^, A^] = and the Hermitian of A^' , so 
Qixv is real. Then, the curvature 2-form can be expressed 
as 

a = -i{A^ A A^) p (16) 

with A® = Afd\». 

A motivation of such separation may be seen from time 
dependence of and if we consider the time evo- 
lution operator U(X) = cxp [— itH(X)] |32j . ft = 1 is as- 
sumed hereafter. We assume the Hamiltonian H(X) and 
the reference state p(X) can be simultaneously diagonal- 
ized by unitary matrix 5(A), i.e., H d (X) = 5 '(A) H (A) 5 f (A) 
and pa — S(X)p(X)S^ (A) are both of diagonal form. 
Correspondingly, we have U(X) = (X)Ud{X)S(X) with 
Ud{X) = exp [— itH d (X)]. Substituting the time evolution 
form of U(X) into Eq. (fT5)l . we obtain 



(i) 



(dpU\d v u) 



(dpU\u)(u\d v u) 



E 

i 



where = ^^{d^V^v^vU^Vv^ and 

(A\B)i = Tr i[piA^B], the subscript I means the l-th sub- 
space. In such cases, we have 



i 



(12) 



with Q\ IV = (dJJ^dJJfo - (fi^l/i|^>i(^|^Ci>l- Hence, 
for the model whose unitary operators and the reference 
state p have the same factorization structure, the OQGT 
is additive. 

Splitting for stationary reference states. For matrix 
which can be diagonalized, the variance with respect to 
the parameters can be decomposed into two parts, the 
variance of the eigenvalues and of eigenvectors. A unitary 
operator U{X) can always be written in the form U(X) = 
S' (A)E/d(A)5(A), where U d (X) is a diagonal unitary matrix 
and 5(A) is unitary. So the derivative of U can be splittcd 
into two terms as follows 



dpU 



(13) 

with = 5 t (a A1 C/ d )5 and D {2) = -i[U,a^], where 

a„ = iS^dnS is a Hermitian matrix. Note that we have 



[d£\u] = and [Di 1] \u] = 0, since D^' and U can 
be simultaneously diagonalized by the unitary operator 

A (2) with 



(i) 



(i) 



5(A). Correspondingly, we have A^ = A r 

Af = itfDf for i = 1, 2. 

For cases in which the reference states p(X) is stationary 
under the unitary operation U (A), i.e., [p, U] = 0, we chose 
5(A) as the unitary matrix diagonalizing simultaneously 
U(X) and p{X). Then we have 

(4 2) ) 



(4 1)f 4 2) ) 



Tr{U*[U,a ll ]p} = 0, 
-zTr {D«t [U, a,] p) 



= 0. 



(14) 



The second equality of second line of the above equations 
is due to the commutation relation [D^\u\ = and 



[p, U] = 0. Then substituting A^ = A^+A^' into Eq. © 
and combining Eq. (|14[) . we obtain the splittcd form of the 
OQGT as follows 

Q$ = (AUAU) P -(AW) P (AW) P , (15) 



{1) 
{2) 



a r 



= E^i 1 - 008 ^-^')*]}, (17) 



where Ei is the eigenvalue of H. a^ v , P 1 ^ are both time- 
independent and defined by 



Tr[(dpH d ){d u H d ) Pd ] 
-TrP M H d )p d ]Tr[(^H d )p d ] 



(18) 



PL = tpM 



where [tt^ij is the matrix element of a M in the represen- 
tation whose basis are the eigenvectors of H. So for such 
cases, Qp 1 } is proportional to i 2 , while Q^ 2 } consists of cir- 
cular functions. It is remarkable that for pure states, a Mt , 
is vanished, only the second term of Eq. (fl~5|) exists. 



XY model. — The one-dimension spin-1/2 XY model, 
with N = 2M + 1 spins, in the presence of a trans- 
verse magnetic field characterized by A, is described by 
the Hamiltonian 



M 



H 



XY 



= - E 



l=-M 



1 



■7 



1 



i V V 



X<7? 



(19) 

where erf (a = x, y, z) are the Pauli matrix for the i-th spin 
and 7 represents the anisotropy in the x — y plane. The 
periodic condition (M + 1 = —M) is assumed here. We 
are interested in a new family of Hamiltonians obtained 
by applying a rotation of <f> around the z direction to each 
spin 

ff(M,A) =ff + (0)ffxY(7,%M, (20) 

where g{4>) = YlfL-M ex P (^4" J iz/2). This model was used 
to establish a relation between geometric phases and crit- 
icality of spin chains [4"0"H4"2~| . 

To obtain the analytic solution of this model, we ap- 
ply three steps of transformations as follows: (i) the 



Jordan- Wigncr transformation 
2c]c 3 )ci, trf = 1 - 2c\a; (ii) 



Fourier transformation 
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q = {l/VN)J2 k e a '" kl/N dk with k = —M, ■ ■ ■ ,M; (iii) 
the pseudo Pauli operator transformation defined by [H] 



nents of the OQGT as follows 



o-fcx = dl<T_ k + d-kdk, 
&ky = —id' k d'_ k + id-kdk, 
o-kz = d\d k + d\ k d- k - 1, 



(21) 



which satisfy the commutative relations {cr ka , <J kb \ — 2<5 a t, 
and [a ka ,cr kb ] = 2ie abc a kc , where e abc is the Levi-Civita 
symbol. For k ^ 0, each of them acts on the Fock space 
spanned by {\v&c}k,- k , \k, —k), \k), \ — k)} with |vac) the 
vacuum state for d k and d- k fermions. These pseudo Pauli 
operators become conventional Pauli matrix in the sub- 
space spanned by {\0} kt _ k , \k, —k}} and vanish elsewhere. 
After these three steps, we get 

M 

H = (\-l)a 0z + H ^ 
k=-M 

H kz = Sl{4>, e k )H k , d (A k )S k (4>, e k ), (22) 

where H kA = A k a kz and S k (4>,6 k ) = R kx {0 k )R kz (4>) with 
R k a(a) = exp(— iaa ka /2) for a = {x, y, z}. The interme- 
diate parameters A k and 9 k are given by 



A* 



l\J (A - cos (2irk/N)) 2 + 7 2 sin 2 {2nk/N), 

^ i A-cos(27rfc/iV) + i7sin(27rfc/iV) 
k ~ ~2 n A - cos (2nk/N) — 17 sin (2nk/N) ' ^ ' 



We consider the OQGT of the time evolution operators 
U = cxp(-itH) and chose the ground state p = \G)(G\ 
as the reference state. They are of the same factorization 
structure 

M 

U = <g) St(<t>, e k ) exp(-itH k!d )S k (cf>, e k ), 

k =0 
M 

P = <g)Sl(6 k )\l) k (l\S k (e k ), (24) 

where | \) k is the eigenstate of a kz with eigenvalue — 1. 
Due to this factorization structure, the OQGT is additive. 
Because [p, U] = and p is pure state, Q^ v is splitted 



(2) 



and Q}j,u of Eq. (|T7)) vanishes. To obtain the Q 
first get the 1-form matrix a = iS^dS, which is given by 
a = Efcio a k with 



a k 



isl{^e k )d[s k ^,e k )] 



Substituting Eq. (|2"5|) into Eq. (fTT)) . we obtain the compo- 



d(j) 



Qxx 

Qsyy 



M 

^ -p- sin 2 (A fe i) sin 2 6 k , 
fe=l k 

M 4 

-p- sin 2 (A fe i) cos 2 6> fc sin 2 
fe=i k 

M 



/2nk\ 



sin 2 (A fc t) [sin(2^)] 2 , 

fe=i 

E4 . 2 . . . / 27rfc\ 

-■p-sm (A fe t) sin cos sin I I 

^ -i— sin 2 (Ajfei) sin 6> fc sin (20 fc ) , 



fe=i 

M 



™ 2 /27rfc 
Q-r^ = X^a - Sin2 ( A fc i ) cos6 'fe sin ( 26 'fc) sm (-jy- 
fc=i fe ^ 

With the help of the metric tensor = ReQ M „, 

we can obtain the Loschmidt echo L = \(G(x)\W (x + 
Sx)U(x)\G(x)}\ 2 for arbitrary multi-parameter perturba- 
tions Sx = (■ • ■ , 5x^ , ■ • • ). This can be seen from Eq. ([2]), 
then we have L = 1 — g^ v 8x^8x v /4. So the metric tensor 
directly reflects the sensitivity of the time evolution to the 
perturbations of multi parameters. A relation between the 
Loschmidt echo and quantum criticality was established in 
Ref. gS]. 

The critical points of the XY model are given by the 
lines A = ±1 and by the segment |A| < 1, 7 = |46j . 
We consider the region of the Ising-typc phase transition 
at A = 1. The time evolutions of the components of the 
metric tensor = Re Q^ v for different values of A are 
shown in Fig. ([T]) , at 7 = 1 for simplicity. The behaviors 
of the A-related components g\\ and gx-y are obviously 
singular at the critical point A = 1. g\\ abruptly increases 
with time near the critical point, which implies the time 
evolution is highly sensitive in this critical point. Near 
the critical point, g\ 7 also has a qualitative distinction, 
it sharply decreases to negative values with time when A 
slightly greater than 1. 

In the three-dimensional manifold of the parameter 
(A, 7, 0), the curvature 2-form a is equivalent to a effective 
vector field B = (B\, B 7 , B^) through the relation 

a = BxdX 1 A dX^ + B 1 dX t ' A d\ x + S dA A A dX 1 . (26) 

In this rotated XY model, we obtain B = 
(2ImQ 7 0, — 2ImQA0, 0). For the case of 7 = 1, the 
time evolution of _B 7 (— 2a\ c f ) ) is very different between 
the two sides of the critical point, as shown in Fig. (J2J. 
When A is slightly less than 1, B 1 sharply increases with 
time, meanwhile when A is slightly greater than 1, £? 7 
sharply decreases from zero to negative values. So we 
conclude that the amplitude of £> 7 sharply increases near 
the critical but the direction of B 1 is opposite in the 
different phases. 
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Fig. 1: Three dimensional diagram of the components of the rescaled metric tensor g^ v — Re Q^/N as functions of A and t for 
the system at 7 = 1 with N = 1001 the number of spins. Here n,v are chosen from the parameter set {A, 7, (fi}. 




Fig. 2: Three dimensional diagram of the components of the rescaled curvature tensor <t m „ = ImQ^/N as functions of A and t 
for the system at 7 = 1 with TV = 1001 the number of spins. Here [i,v are chosen from the parameter set {A, 7, <fi}. 
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Conclusion. — To summarize, we have introduced 
the operator-quantum-geometric tensor on the manifold 
of the unitary operators with a given reference state, and 
prove the additivity for factorizable models. We show the 
splitting of two kinds contributions to the OQGT for the 
cases of stationary reference states, and these two con- 
tributions have different type of time dependence. This 
OQGT approach can be applied to investigate the sen- 
sitivity of unitary operation against perturbations. The 
Loschmidt echo for the unperturbed and perturbed time 
evolutions for arbitrary kind of multi-parameter perturba- 
tions can be easily got as long as the OQGT is obtained. 
We used this approach to investigate sensitivity of the time 
evolutions of the rotated XY models and show the sharp 
changes of the OQGT near the quantum critical points. 
We believe that this OQGT approach is useful to some 
interesting questions, like the quantum criticality pTl,33 , 
decoherence [35] and quantum chaos [MJI35], e ^ c - After 
the acceptance of this paper, we notice the work |47j where 
both the QGT and Operator fidelity are shown to be rele- 
vant to the adiabatic error in the adiabatic and holonomic 
quantum computing. 
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